What controls the Fe ii strength in active galactic nuclei? 
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ABSTRACT 

We have investigated the correlations of the equivalent widths (EWs) of narrow and 
broad Fe II emission lines and the Fe II/Mg II intensity ratio with fundamental physical 
parameters of active galactic nuclei (AGNs), luminosity (L), black hole mass (Mbh) 
and Eddington ratio (L/LEdd), using a homogenous sample of z < 0.8 Seyfert 1 galax- 
ies and QSOs in the spectroscopic data set of Sloan Digital Sky Survey Fourth Data 
Release (SDSS DR4). The sample comprises 4178 spectra that are selected to suffer 
little from the host-galaxy starlight contamination. We find that the strongest correla- 
tions of almost all the emission-line intensity ratios and EWs are with L/L^m, either 
positively (e.g. Fell EW) or negatively (e.g. Mgll EW), rather than with L or Mbh; 
In particular, the intensity ratios of both the ultraviolet and optical Fe II emissions to 
Mgll A2800 correlate quite tightly with L/Lndd (e.g.. Spearman = 0.74 for the [nar- 
row FeII]/MgII - L/Ledd correlation). We have also investigated the Fell emission in 
the full sample of about 27,000 type 2 AGNs in the SDSS DR4, and find that no Fell 
is detected at 3a significance. We argue that (1) the narrow-line Fell in AGNs origi- 
nates purely from the inner narrow-line region on scales less than the torus whereas the 
sites of the broad-line Fe II are likely not unique; (2) The variation of the emission-line 
strength in AGNs are regulated by L/L^dd because it governs the global distributions 
of the properties (particularly column density) of the clouds gravitationally bound in 
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the line-emitting region; (3) The systematic dependence on L/L^^ must be corrected 
for when using the Fe Il/Mg II intensity ratio as a measure of the Fe/Mg abundance 
ratio to study the history of chemical evolution of quasar environment. 

Subject headings: radiation, accretion - galaxies: active - quasars: emission lines - 
quasars: general 

1. Introduction 

Fell multiplet emission is a prominent feature in the ultraviolet (UV) to optical spectra of 
most active galactic nuclei (AGNs), including QSOs at redshifts as high as > 6 (e.g., Freudling et 
al. 2003, Barth et al. 2003, Iwamuro et al. 2004, Jiang et al. 2007, Kurk et al. 2007). Its equivalent 
width (EW) varies significantly from object to object (e.g., from > 100 A to undetectable for the 
emission in the 4434^684A region [hereafter A4570]; Boroson & Green 1992, hereafter BG92). 
Fell is an important probe of AGN physics. For example, BG92 showed that the relative strength 
of the Fe II A4570 emission is a dominating variable in the first principal component (PCI) in their 
analysis of the correlation matrix of QSO properties which is generally believed to be linked to cer- 
tain fundamental parameter of accretion process (e.g., the relative accretion rate, often expressed 
as L/L^dd or L/M; Sulentic et al. 2000, Boroson 2002). More importantly. Fell emission relative 
to that of a-elements (e.g. Mg II) can be considered to be an observational proxy of the Fe/a abun- 
dance ratio. If iron is mainly produced by Type la supernovae that have long-lived progenitors, 
this ratio would have a characteristic delay of ~ 1 Gyr from the initial starburst in the quasar host 
environment, based on chemical evolution models. Thus Fell emission in principle can serve as a 
cosmic clock to constrain the age of QSOs and the epoch of the first star formation of their host 
galaxies (Hamann & Ferland 1993; and cf., e.g., Matteucci & Recchi 2001). 

However, the ionization and excitation mechanisms of the Fe II emission are very complex due 
to the low ionization potential of Fe° atom (7.9 eV), the low energy levels of various excited sates 
of Fe^ ion, and particularly the complexity of its atomic structure (see Baldwin et al. 2004, Collin 
& Joly 2000 and references therein, Joly et al. 2008). Thus the conversion from Fell emission to 
iron abundance is not straightforward (e.g. Verner et al. 2004, Baldwin et al. 2004). Furthermore, 
it is possible that Fell emission of AGNs arises from a number of different sites with similarly 
suitable excitation conditions, including gaseous clouds gravitationally bound in the narrow-line 
and broad-line region (NLR and BLR), (the base of) outflows, and the surface of the accretion disk 
(AD) [see, e.g., CoUin-Souffrin 1987 et seq., Murray & Chiang 1997, Matsuoka et al. 2007, Zhang 
et al. 2007]. This has already been hinted by the reverberation observations of Fell emission, 
particularly of its optical component, since the cross-correlation function appears to be flat-topped 
(Vestergarrd & Peterson 2005, Kuehn et al. 2008; see also KoUatschny & Welsh 2001, Wang et 
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al. 2005). Therefore, to use of Fe II emission as a proxy of Fe abundance and further as a cosmic 
clock requires a better understanding of its emission mechanism. In fact, it is well known that the 
Fe II/Mg II ratios of QSOs at the same redshift have a rather large scatter (e.g. Iwamuro et al. 2004, 
Leighly & Moore 2006). 

It is also possible that the excitation mechanisms and sites of Fell emission are well gov- 
erned by certain fundamental parameters (such as Eddington ratio £] due to certain self -regulation 
mechanisms that maintain the normal dynamically quasi-steady (or even quasi-virialized) states of 
the gas surrounding the central engine of AGNs. It is well-known for decades that QSO spectra 
are remarkably similar, at the zeroth-order approximation, with emission lines of similar inten- 
sity ratios sitting atop a blue continuum regardless of their luminosities and redshifts (Davidson 
& Netzer 1979, Korista 1999, Dong et al. 2008; see more evidence and insightful arguments for 
the similarity in Gaskell et al. 2004, Gaskell et al. 2007, Gaskell & Benker 2008)0 This sim- 
ilarity means that the integrated spectra and global properties of QSOs are rather insensitive to 
the detailed distributions of various gas properties and the 'microscopic' physical processes taking 
placing in the line-emitting region of AGNs; instead, they are governed by the powerful selection 
mechanisms introduced by the atomic physics (albeit not exactly identified yet), as proposed by 
Baldwin et al. (1995) in the "locally optimally-emitting clouds" model that has reproduced the 
overall QSO broad-line spectra (e.g. Korista & Goad 2000). If the variations of the emission lines, 
particularly, the narrow-line or broad-line Fell emissions, are in a similar situation — i.e., their 
variations are indeed governed by certain self -regulation mechanisms and thus are regular, we may 
calibrate the translation from Fe II emission to abundance and constrain chemical evolution in the 
quasar environment. 

Motivated by these considerations, we carry out a systematic study on the strengths of UV 



' Eddington ratio {£ = L/LEdd) is the ratio between the bolometric and Eddington luminosities. Eddington lu- 
minosity (LEdd), by definition, is the luminosity at which the gravity of the central source acting on an electron- 
proton pair (i.e. fully ionized gas) is balanced by the radiation pressure due to electron Thomson scattering; 
^Edd = 47rGcMmp/o'T, where G, c, M, nip , crj are the gravitational constant, speed of light, mass of the central source, 
proton mass, Thomson scattering cross-section, respectively. In accretion-powered radiation systems, L/^Edd is often 
referred to as dimensionless accretion rate m (the relative accretion rate normalized by Eddington accretion rate MEdd, 
rh = M/MEdd = / L^ii, M being mass accretion rate and 77 the accretion efficiency) as m is not an observable; yet 
the two notations are different both in meaning and in scope of application. Even in the accretion-powered radiation 
systems like AGNs, L/^Edd (X) is not equivalent to m (M) except in the simple thin accretion disk model of Shakura & 
Sunyaev (1973). Therefore, we call L/LEdd dimensionless luminosity (£). For clarification or for the ease of notation, 
below we use L/LEdd and £ alternately. 

^The zeroth-order similarity perhaps holds even for low-luminosity and low-accretion AGNs (see, e.g., Moran et 
al. 1999 for NGC4395, Maoz 2007). If it is the case, it means that the global spectral energy distribution and emission 
Unes of AGNs is quite insensitive to Mbh in the range of 10^ - 10'" Mq and L/L^dd in the range of 0.001 - 1 . 
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and optical Fe II emissions, by taking advantage of the unprecedented spectroscopic data from the 
Sloan Digital Sky Survey (York et al. 2000). Particularly, using the analysis of Veron-Cetty et 
al. (2004), we are able to study narrow-line Fell emission systematically for the first time. The 
paper is organized as follows: In section 2 we describe the selection criteria of the sample, the 
spectral processing, and correlation and regression analyses. Section 3 presents the results and our 
discussions. Section 4 is our brief conclusions and implications. Throughout this paper, we use a 
cosmology with //o=70 kms"' Mpc"\ fiM=0.3 and Qa=0.7. 



2. Data analysis 
2.1. Sample and measurements 

2.1.1. Sample construction 

We first construct a homogenous sample of Seyfert Is and QSOs (namely type-1 AGNs) from 
the spectral data set of the SDSS Fourth Data Release (Adelman-McCarthy et al. 2006), according 
to the following criteria: (1) redshift z < 0.8; (2) the median signal-to-noise ratio (S/N) > 10 per 
pixel in the optical Fe ll and H/3 region (4400-5400A); (3) the rest-frame absorption-line equivalent 
widths (EWs) of Ca K (3934A), Ca K -i- He (3970A) and (4102A) < la. The redshift cut is set to 
ensure the H/? and Fe II A4570 multiplets are present in the SDSS spectra; the S/N criterion allows 
proper placement of local continua and the accurate measurement of emission lines (especially 
broad-line and narrow-line Fell emissions); the third criterion ensures the measurement of the 
AGN luminosity and the emission-line EWs suffering little from the contamination of the host- 
galaxy starlight. By combining SDSS spectra of high-luminosity QSOs and absorption-line or 
star- forming galaxies we find that our selection criteria corresponds to a galaxy contribution < 10% 
around 5 100 A. 

After removing duplications and sources with large number of bad pixels in the H/? and Fe II 

A4570 region, our final sample consists of 4178 type-1 AGNs (hereafter the full type 1 sample). 
Of them, 2092 have redshift z > 0.45 and median S/N > 10 per pixel in the UV Fe II region (2200- 
3090A), which are used to examine the behavior of UV Fe II and Mg II AA2796, 2803 doublet lines 
(hereafter A2800), and referred to as the UV subsample. 



2.1.2. Spectral fitting and measurements 

We use the measured parameters from the Valued-added ExtraGAlactic Catalog developed 
and maintained by Center for Astrophysics, University of Science and Technology of China (USTC- 
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VEGAC; X.B. Dong et al., in preparation, J.G. Wang et al., in preparation). Details of spectral anal- 
ysis will be described in separate papers. We discuss briefly here the steps relevant to this paper: 
we correct the SDSS spectra for Galactic extinction, separate Fe II emission from the continuum, 
and fit the emission lines. The optical Fell is modeled with two separate sets of templates in an- 
alytical forms, one for the broad-line system and the other for narrow-line system, constructed 
from the measurements of IZW 1 by Veron-Cetty et al. (2004) [see Dong et al. 2008 for details]. 
The regions of UV broad Fe n is modeled with the tabulated semiempirical template generated by 
Tsuzuki et al. (2006) basically from the measurements of I ZW 1 ; in the wavelength region covered 
by Mg II emission, this template also employs the calculation with the CLOUDY photoionization 
code (Ferland et al. 1998). Mgll doublet lines each is fitted with a truncated 5-parameter Gauss- 
Hermite series (Wang et al. 2009; see also Salvainder et al. 2007). As its profile is sometimes 
rather complex, broad-line H/5 is fitted with as many Gaussians as statistically justified; we do not 
ascribe any particular physical significance to the individual Gaussian components (Dong et al. 
2005). Every broad Fell lines are assumed to have the same profile as broad H/9 while the redshift 
of broad-line Fe II is set to be a free parameter; the uncertainty of the Fe II line width has little 
effect on the measurement of the Fell flux (Salvainder et al. 2007, Hu et al. 2008). All narrow 
emission lines are fitted with a single Gaussian except the [O III] A A4959,5007 doublet; each line 
of the [O III] doublet is modeled with two Gaussians, with one accounting for the line core and the 
other for a possible blue wing as seen in many objects. The redshift and width of narrow-line Fe II 
are set to be free parameters. 

Here we must note that, to properly measure the parameters of Fe II emission, the separation 
of narrow-line Fe II from the broad-line one is very important; as far as we know, this point has 
been long neglected mainly because of the lack of a narrow-line Fe II template until Veron-Cetty et 
al. (2004). Narrow Fell multiplets, both permitted and forbidden, have been observed in the AGN 
spectra, both in the UV (e.g., Vestergaard & Wilkes 2001, Wang et al. 2008) and in the optical 
(e.g., Veron-Cetty et al. 2004, Veron-Cetty et al. 2007, Bruhweiler & Verner 2008; see also §4.6 
of Hu et al. 2008). Certainly, whether there is a physically distinct boundary between the NLR and 
BLR of AGNs is still hotly debated (cf. Laor 2007), particularly considering the likely presence 
of the region intermediate between them (e.g., the Intermediate-Line Region [ILR], Brotherton et 
al. 1994; the inner NLR, Nagao et al. 2003, see Zhang et al. 2009 for a concise review); and 
both the BLR and NLR may be highly stratified and multi-zoned. In this paper, we follow Veron- 
Cetty et al. (2004) who phenomenologically took the Fe II component presented in other broad 
lines as the broad-line (LI), the Fe II (both permitted and forbidden) component presented in other 
low-excitation narrow lines the narrow-line (N3). In IZW 1, there are two additional systems, Nl 
and N2, that appears only in high-excitation lines (not in Fell); Nl and N2 are relatively broad 
and blueshifted by 1450 and 500 km s~^ respectively whereas N3 is almost at the systemic redshift 
of the host galaxy. The broad-line Fell system of IZW 1 is also blueshifted by about 150 - 200 
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km in both the optical and UV bands (Veron-Cetty et al. 2004, Laor et al. 1997; and Cf. Hu 
et al. 2008), which distinguishes it from the narrow-line system. The in-depth investigation on the 
physics of the narrow-line Fell is beyond the scope of this paper and deferred to a forthcoming 
paper. 

In Figure 1, we show the optical spectrum of a representative object in the present sample 
(SDSSJ092801. 30+491817.3), as well as the result of our spectral decomposition. In this object, 
the optical Fell broad lines have a width of 1300 km s"' FWHM; narrow Fell, 250 km s"' (cor- 
rected for the instrumental broadening of 130 km s"' FWHM). Note that the individual narrow 
Fe II lines marked in the figure (particularly the forbidden lines) are separate from any broad Fe II 
lines, which indicates that our decomposed narrow Fell multiplets are bona fide, not residuals 
from broad Fe II subtraction due to mismatch of the broad-line Fe II model. The AGNs with strong 
narrow-line Fell emission are of diverse types: NLSls and non-NLSls, broad absorption-line 
(BAL) QSOs and non-BAL QSOs, radio-loud and radio-quiet. Figure 2 displays several diverse 
SDSS spectra with strong narrow Fe II emission in the sample, as well as our decomposition in the 
rest-wavelength range of 4000-5600A; in its bottom panel, a low-ionization BAL QSO with vari- 
able, high radio luminosity (FBQS J152350.4+391405, Ljghz = 4.0 x 10^^ ergs s"^ Hz-^; Becker et 
al. 2000) is displayed. In Figure 3, we show the spectral fitting in the UV Fell region. 

For all measured emission-line fluxes, we regard the value as reliable when it is > 2a; other- 
wise, we adopt the value plus the 2a error as an upper limit. We calculate the black hole masses 
based on H/? line using the relation presented in Vestergaard & Peterson (2006, their Eqn. 5). 
This relation was calibrated with reverberation mapping-based masses and used the R-L relation 
by Bentz et al. (2006) which was corrected for host-galaxy starlight contamination. The Edding- 
ton ratios are calculated assuming a bolometric luminosity correction Lboi ~ 9ALa(5100) (Kaspi 
et al. 2000, Elvis et al. 1994). The mean and standard deviation (computed in log-space) of the 
FWHM of broad H/? is 3600 km s^^ and 0.22 dex; ALa(5100), 4 x lO^^ergs s'^ and 0.40 dex; Mbh, 
2 X IO^Mq and 0.52 dex; L/L^dd, 0.13 and 0.44 dex. Because the SDSS spectroscopical survey is 
magnitude-limited, Mbh and L/Lndd correlate with each other apparently; this apparent correlation 
is enhanced by the mutual dependence of both variables, by construction, on FWHM of broad H/? 
and ALa(5100). The Spearman correlation coefficients (r^) between Mbh and L/L^dd are -0.663 
for the full type 1 sample (4178 objects), and -0.762 for the UV subsample (2092 objects). 



2.2. Correlation and regression 

We use the generalized Spearman rank method implemented in the ASURV package (Isobe et 
al. 1986) to perform bivariate correlation test, and the generalized partial Spearman rank method to 
perform partial correlation test (Kendall & Stuart 1979; see, e.g., Shapley et al. 2001 for a previous 
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application). Both methods test for not only a linear relation but a monotonic one; and they can 
handle censored data in both independent and dependent variables. The results are summarized in 
Table 1, where we report the Spearman coefficient (r^) and the probability (Pnuii) that a correlation 
is not present. Out of the correlation results, two striking features emerge as follows: 

• According to the Spearman correlation analysis, the strongest correlations of almost all 
the emission-line intensity ratios and EWs are with L/L^dd, rather than with L or Mbh it- 
self! This is confirmed by the Spearman partial correlation analysis: when controlling for 
linewidth, L or Mbh, there are still strong correlations with L/L^dd, but NOT vice versa. Be- 
sides, generally line ratios have tighter correlations with L/L^dd than the EWs of the related 
emission lines. 

• The intensity ratios of Fe II emissions - both narrow and broad, both in the UV and in the 
optical - to Mgll A2800 have very strong, positive correlations with L/L^dd- Interestingly, 
(narrow FeII)/MgII has a stronger correlation with L/L^dd than (broad FeII)/MgII (r^= 0.74 
versus 0.58 and 0.46). 

In Figure 4, we plot the equivalent widths of the ultraviolet (2200-3090A), optical (4434- 
4684A) broad-line and narrow-line Fell emissions versus ALa(5100A), black hole mass, and Ed- 
dington ratio {i = L^^oi/L^dd)- In Figure 5, the three correlations of Fell/Mgll ratios with L/L^dd 
are illustrated; these tight correlations are striking considering the narrow range of L/L^dd in our 
sample (Icr = 0.44 for a log-normal distribution) and the known systematical errors in the L/L^dd 
estimation (see §3.1.1). We performed linear regressions in the log-log scale using the method 
of Kelly (2007). This method adopts a Bayesian approach and accounts for measurement errors, 
censoring and intrinsic scatter. The results are as follows, 

(-0.18 ±0.05) + (1.77 ±0.05)log^ (1) 

(0.34 ± 0.02) + (0.74 ±0.03) log ^ (2) 

(0.94 ± 0.01) + (0.33 ±0.02) log£ (3) 

^The only exception is the EW of broad H/3. This is because a significant fraction of broad H/3 is emitted from the 
surface of the AD (and/or high-ionization, optically-thin clouds). We will discuss this point and other properties of 
broad H/3 in a forthcoming paper (see also §3.3). Note that Hydrogen Balmer lines, which can be produced through 
both the recombination process of ion and the coUisional-excitation of H" atom, originates efficiently from both 
the fully ionized H""" region and the partially ionized HI* region; thus the total emission of Balmer lines is hybrid, 
as suggested by mounting observational facts (see Marziani et al. 2008 and references therein; §4.2 of Wang et al. 
[2009]). 
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The intrinsic standard deviation of these relations (corrected for measurement errors) are 0.04, 0.22 
and 0.15, respectively. 



3. Results and Discussions 
3.1. L/LEdd controls Fen strength 

3.1.1. The phenomenological regularities 

As shown above, a general trend is that the emission-line EWs and intensity ratios correlate 
more tightly with L/Lsdd than with L or Mbh- In this paper we focus on the behaviors of the 
narrow-line and broad-line Fe n emissions, and defer the comprehensive discussions to an accom- 
panying paper (Dong et al., 2009b). First, it is quite unexpected that narrow-line Fell has a tighter 
correlation with L/Lsdd than broad-line Fell, as the gas emitting broad-line Fell is thought to be 
closer to and thus should be linked more tightly to the central engine than the narrow-line. One 
possible explanation is that the origin of narrow-line Fe II is more homogeneous than the broad- 
line (see §3.2 and 3.3). Second, the line ratios have stronger correlations with L/L^dd than the 
EWs of narrow and broad Fe II. This is natural since, compared to line ratio, EW depends on more 
additional parameters such as the covering factor of the BLR and the anisotropy of the continuum 
radiation. The exceptions (see Table 1) are the correlations between line ratios and L/L^d that 
have a flat log-log slope, such as (narrow Fell)/(broad Fell); the dynamical range of the line ratio 
is small and the relation is dominated by fluctuation caused by the measurement errors. Third, 
UV Fe II has not as tight a correlation with L/Lsdd as the optical Fe ll. This is caused by certain 
radiation transfer effects that we will discussed in Dong et al. (2009b). Fourth and probably the 
most important, as we mentioned in §2.2, in any case the ratios of Fell to Mgll have quite tight 
correlations with L/L^dd, and yet not with L or Mbh- 

Collin et al. (2006) argued that the use of FWHM would introduce systematic biases to 
the virial Mbh estimate due to viewing angle effect and the presence of non-virialized velocity 
component(s) (see also Sulentic et al. 2006; Wang et al. 2009). We thus calculate Mbh also using 
the relation proposed by Collin et al. (2006, their Eqn. 7) and others; this alternative Mbh estimate 
leads to similar results. We note that, in the plot of (narrow Fell)/Mgll versus L/L^dd (panel c of 
Figure 4), there are some data points at the low L/L^dd end (< 0.05) with systematically larger line 
ratios. This is likely caused by the over-estimated Mbh due to the systematic biases as suggested 
by those authors, which is confirmed by the correlation of the offset with FWHM (r^ = 0.46 for the 
1492 objects with reliably detected narrow Fe II). All the outliers with offsets of > 1 dex upward 
have FWHM > 5500 km s"' (with the mean and standard deviation of 9200 and 2200 km s~^). 
For such large FWHMs, even Collin et al.'s correction might be inadequate (in fact, Collin et 
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al.'s correction is by a factor of only 0.5 compared to the estimates of the Vestergaard & Peterson 
formalism). We will address this issue using a well-calibrated Mbh formalism with Mgll FWHM 
in Wang et al. (2009). 



3.1.2. The physical picture 

A strong, negative correlation between EW and L/L^m has been noted by several authors 
for Civ A1549 line (Baskin & Laor 2004; see also Bachev et al. 2004). Both groups suggested 
that L/Ledd is the fundamental driver of the Baldwin effect (BEff) - a well-known anticorrelation 
between the emission-line EWs and the AGN luminosity (Baldwin 1977), although they then did 
not know the underlying physical mechanism. The findings in this paper expand this picture. 
Although the gas environment in the line-emitting region of AGNs is complex and chaotic, the 
strength of almost all emission lines are governed by L/Ledd- It is well-known that C IV is a high- 
ionization line; Mgll, low-ionization and optically thick; and Fell, low-ionization and optically 
thin (e.g., CoUin-Souffrin et al. 1986). The emerging pattern is that, as L/L^dd increases, the 
emission strength of high-ionization or optically thick lines decreases whereas the strength of low- 
ionization and optically thin lines increases. 

High-ionization lines are emitted from the illuminated surface of the clouds; Optically thick 
lines, either from the illuminated surface (e.g. recombination line Lya) or from the thin transition 
layer located immediately behind the Hydrogen ionization front (e.g., with the optical depth at the 
Lyman limit being in the range of 10 and 10"^ for low-ionization Mg II; CoUin-Souffrin et al. 1986); 
yet low-ionization, optically thin lines such as Fe II multiplets arise from the partially ionized H I* 
region (i.e., from the ionization-bounded clouds only). Hence, the correlation pattern may be 
telling us that, as L/L^dd increases, the column density (A'h) of the clouds in the line-emitting region 
increases. We argue that the underlying physics is certain self-regulation mechanisms caused by (or 
corresponding to) Eddington ratio (dimensionless luminosity, not dimensionless accretion rate)0 
These mechanisms maintain the normal dynamically quasi-steady (or even quasi-virialized) states 
of the gas surrounding the central engine of AGNs; i.e., they regulate the global distributions of the 
properties of the clouds gravitationally bound in the line-emitting region (including the inner NLR, 
see §3.2). The details are beyond the scope of this paper, and are deferred to Dong et al. (2009b). 



Certainly, for particular radiation systems that are powered by gravitational accretion (e.g. AGNs), £ is linked 
tightly to m anyway (cf. Footnote 1; Merloni & Heinz 2008). Thus there is an additional effect associated with £ yet 
directly via m as follows. The increase in £ means, meanwhile, the increase in m, the gas supply. Reasonably, it is from 
the supplied gas spiraling into the central engine that (at least a significant part of) the line-emitting clouds originate; 
this is particularly true for the BLR and inner NLR clouds that are located between the torus (as the fuel reservoir) and 
the accretion disk (e.g. Gaskell & Goosman 2008). Hence, as £ increases, the total mass of line-emitting gas increases. 
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Essentially, there is a lower limit on the column density of clouds bound to the AGN line-emitting 
region, which is set by L/L^dd- Even at L/L^dd < 1, Iow-A^h clouds are blown away because they 
are not massive enough to balance the radiation pressure force that is boosted by photoelectric 
absorption (see, e.g., Fabian et al. 2006, Marconi et al. 2008). According to Figure 1 of Fabian et 
al. (2006), for dust-free clouds of A^h ^ 10^^ cm"^, the low limit is approximately proportional to 
L/LEdd as Nu > 10^^ i cm"^; for dusty clouds, A^h > 5 x 10^^ i cm"^. 

Last, we note that the tight, positive Fell/Mgll - L/L^dd correlations cannot be caused (or 
enhanced) by any intrinsic relation between the Fe/Mg abundance ratio and L/L^dd due to the 
possible starburst-AGN connection. Because in any plausible scenario of the connection that the 
AGN activity is triggered by the starburst (e.g. Sanders et al. 1988, Davies et al. 2007), no matter 
how long is the delay between the two events, the AGN is first fully accreting and then gradually 
weakening; this would imply that Fe/Mg abundance ratio would correlates negatively with L/L^dd 
opposite to the trend we see here. 

3.2. The site of the narrow-line Fe n emission 

As discussed above, the tight correlation between narrow-line Fell and L/L^dd hints that the 
gas emitting narrow-line Fe 11 probably is very homogeneous and thus the emission is not affected 
significantly by additional parameters. To investigate this, we search the Fe II emission in all the 
~ 27,000 type 2 AGNs in the SDSS DR4, which are selected according to the selection criteria of 
Kauffmann et al. (2003) and have S/N > 5 for H^, [Oin] A5007, Ha and [Nil] A6583 lines. The 
mean and standard deviation of their redshifts are 0.10 and 0.05 respectively. No Fe II emission is 
detected at > 3cr significance in any object in the sample. Thus, according to the well-established 
unification scheme for Seyfert galaxies (e.g. Antonucci 1993), we conclude that narrow-line Fell 
originates purely from the gas bound in the inner NLR that locates closer to the central engine 
than the torus. The inner edge of the dusty torus is known to be on scales of parsecs, roughly the 
dust sublimation radius; its extent is likely on scales of several tens of parsecs (e.g., Schmitt et al. 
2003, Jaffe et al. 2004, Zhang et al. 2008). 



3.3. The sites of the broad-line Fe n emission 

In contrast to the narrow-line Fe II, the EWs of the broad-line Fe II emissions show at most 
a moderate correlation with L/L^dd- What caused the different behavior? A natural explanation 
is that broad Fen does not originate from a unique site, i.e., not only from the clouds bound in 
the BLR. An additional site is plausibly the surface of the AD, which has been speculated and 
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argued for a long time (e.g. CoUin-Souffrin et al. 1980, Murray & Chiang 1997, Zhang et al. 
2006). The AD surface at radii of typically hundreds of Tq (ro = GM/c^) is argued to be an ideal 
site to emit low-ionization lines such as Balmer lines, Mgll A2800 and the UV and optical Fell 
(e.g. CoUin-Souffrin 1987 et seq.). This is convincingly supported by the observations of low- 
ionization lines of double-peaked profile in a minority of AGNs (Balmer lines: e.g. Chen et al. 
1989, Eracleous & Halpern 2003, Strateva et al. 2003; Mgll: e.g. Halpem et al. 1996), particularly 
by the discoveries of Balmer lines having, in addition to a symmetric broad component located at 
the system redshift, an separate, extremely broad double-peaked component that is apparently 
gravitationally redshifted (Chen et al. 1989, Wang et al. 2005; see a treatment for the general case 
of twisted, warped disk in Wu et al. 2008). This disk component constitutes partly the often-called 
"red shelf" or "very broad component" (see, e.g., Sulentic et al. 2000) of the Balmer-line profile. 
In fact, double-peaked Fe II emission lines are observed in some disk- accreting systems such as 
cataclysmic variable stars, together with other double-peaked low-ionization lines such as Balmer 
lines, Call, Si II and He I lines that arise definitely from the AD according to Doppler tomography 
(e.g. Roelofs et al. 2006, Smith et al. 2006). 

Note that the AD radii that can emit low-ionization lines are less than the critical radius of the 
gravitationally unstable part of the AD, r^it 2 x ^^(Mbh/IO^Mo)-"-^^^ (Collin & Hure 2001). 
It is argued that the gravitationally unstable disk is one source of the clouds bound in the line- 
emitting region of AGNs, particularly for low-ionization lines such as Mgll and both narrow and 
broad Fe II as suggested in this paper (Collin & Hure 2001 ; see also Leighly 2004); other sources 
of the clouds can be outflows that are liable to fragmentation (Proga et al. 2008) and gas infalling 
from the torus (Gaskell & Goosmann 2008). Compared to the gas in the torus, these clouds have 
less or even no dust. 

Finally, we must mention that, in such high-density and high-column density clouds, some 
researchers have argued that photoionization heating might not be sufficient, additional heating 
like a mechanical one seems be necessary to strengthen the Hi* region (e.g., Veron-Cetty et al. 
2006, Joly et al. 2008). There is evidence that mechanical heating from outflows may be such a 
source (S.-H. Zhang, in preparation), although its contribution is generally expected not to be very 
significant and hard to estimate observationally. Outflows, which are launched from the AD — 
dependent on L/L^dd in a certain way, can have large inclination angles or even equatorial (Murray 
et al. 1995, Proga et al. 2008), and thus collide with clouds in the line-emitting region or even in 
the torus. The collisions would also increase clouds in the BLR and inner NLR. 
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4. Conclusions and implications 

In this paper, we find that the narrow-line Fe II originates purely from the inner NLR on scales 
less than the torus. Yet, the sites of the broad-line Fell are likely to be not unique: similar to the 
situation of broad H/3 (probably to a less degree), a significant fraction of the emission is from the 
surface of the accretion disk instead of the clouds bound in the BLR. In summary, Fe II emission 
will arise from any gas surrounding the central engine of AGNs that has sufficiently high density, 
column density and heating energy input, and at the same time has very low dust content. 

More importantly, although the excitation mechanisms of Fe II emission (and the sites for the 
broad-line one) are complex, its relative strengths with respect to both the continuum and other 
emission lines (particularly Mgll) correlate tightly and positively with L/Ledd, and intrinsically 
not with AGN luminosity or Mbh- That is, in addition to the global similarity (to the Oth-order 
approximation) of the QSO spectra, the Ist-order variation of the emission lines are controlled by 
L/Lndd- The essential underlying physical mechanism of these correlations is that Eddington ratio 
(dimensionless luminosity, i. = L/L^dd, see Footnote 1) regulates the distribution of the column 
density of the clouds bound in the line-emitting region. We will discuss it further in an accompa- 
nying paper, including its implication to the BEff and its second-order effects (Dong et al. 2009b). 
One possible interesting implication is that, as far as spectral correlations in the UV-optical band 
are concerned (particularly the Fell - [Olll] anticorrelation), BG92's PCI could share the same 
origin - L/Lgdd - with PC2, which is exactly the BEff of Hell A4686 line (BG92, Wills et al. 
1999; see Dietrich et al. 2002 for an investigation of the He II BEff). 

If the observed large scatter of Fe Il/Mg II at the same redshift is caused predominately by the 
diversity of L/Ledd, then it is hopeful to still use Fell/Mgll as a measure of the Fe/Mg abundance 
ratio to study the chemical evolutionary history of the gas surrounding AGNs and to constrain its 

age ("cosmic clock"), once the systematic variation caused by L/L-^dd is corrected according to the 
formalisms (Eqn. 1-3) presented in this paper. We will address this issue in a forthcoming paper 
(H.Y. Wang et al., in preparation). 
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Fig. 1. — Demonstration of spectral fitting in the optical region. In the figure, we show: the SDSS 
spectrum (black), the local AGN continuum (blue), the continuum plus narrow-line Fell complex 
emission (red), the continuum plus broad-line Fell complex emission (green), the pseudocontin- 
uum (continuum plus total Fell, navy), and the pseudocontinuum subtracted residual (gray). Also 
marked are some narrow Fell lines (particularly the forbidden lines) that are separate from any 
broad Fell lines. The strong emission lines in the residual spectrum are chopped for clarity. 
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Fig. 2. — Demonstration of the SDSS spectra of several strong narrow Fe U emitting AGNs of 
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Fig. 3. — Demonstration of spectral fitting in the ultraviolet region. In the figure, we show: The 
SDSS spectrum (black), the local AGN continuum (blue), the pseudocontinuum (continuum plus 
Fell complex emission, red), and the pseudocontinuum subtracted residual (gray). 
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Fig. 4. — The relation between the equivalent widths of the ultraviolet (2200-3090A), the optical 
(4434-4684A) broad-line and narrow-line Fen emissions of type-1 AGNs and their ALa(5100A), 
black hole mass (Mbh), and Eddington ratio (i = Lboi/LEdd)- 
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Fig. 5. — The relation between the intensity ratios and Eddington ratio for type-1 AGNs. Also 
plotted are the respective best-fitted linear relation in the log-log scale (see Eqn. 1-3). Note that in 
panel (c) there are some outliers at the low L/Ledd end with systematically larger line ratios (e.g., 
by 1 dex); this is likely caused by the over-estimated Mbh due to some systematic biases because 
all the outliers have FWHM > 5500 km s"' (with the mean and standard deviation being 9200 and 
2200 km s"^). 



Table 1: Results of Spearman correlation and partial correlation analysis ^ 



X 



FWHM (H/3B) 



ALa(5100) 



(X , LI Leu ; Mbh) " (X , Mbh ; i/^Edd) 



EW (Fell-'"' 4570) 
(Fell" 4570)/MgII 
(FeH'^ 4570)/H/3B 
(FenJ^4570)/[O III] 5007 
EW (Fen" 4570) 
(Fen*4570)/Mgn 
(Fell'*4570)/H/3B 
(Fell^ 4570)/[O III] 5007 
EW(FellUV) 
(FeIIUV)/MgII 
(FenUV)/H/3'' 
(Fen UV)/[0 III] 5007 
(Fen^v 4570)/(FeII UV) 
(Fell'^ 4570)/(FeII UV) 
(Fell" 4570)/(Fen* 4570) 
EW([OIII] 5007) 
EW (H/3I*) 
EW(Mgn) 



-0.659 (<le-14) 
-0.731 ( <le-14) 
-0.680 (<le- 14) 
-0.606 (<le-14) 
-0.333 (<le-14) 
-0.567 ( <le-14) 
-0.474 (<le-14) 
-0.244 ( <le-14) 
0.023 ( 3e-01 ) 
-0.429 (<le-14) 
-0.159 (<le-14) 
-0.079 ( 3e-04) 
-0.698 ( <le-14) 
-0.435 ( <le-14) 
-0.584 (<le-14) 
0.104 ( <le-14) 
0.333 ( <le-14) 
0.496 (<le-14) 



0.102 ( <le-14) 
0.097 (<le- 14) 
0.079 (<le- 14) 
0.140 (<le-14) 
0.156 (<le-14) 
0.116 ( <le-14) 
0.064( <le-14) 
0.174 ( <le-14) 
-0.036 ( le-01 ) 
0.157 (<le-14) 
0.044 ( 4e-02 ) 
0.103 (<le- 14) 
0.060 ( 6e-03 ) 
0.042 ( 5e-02) 
0.066 ( 2e-05 ) 
-0.098 (<le- 14) 
0.110 (<le-14) 
-0.216 (<le-14) 



-0.508 (<le-14) 
-0.627 (<le- 14) 
-0.537 (<le-14) 
-0.440 (<le-14) 
-0.230 (<le-14) 
-0.480 (<le- 14) 
-0.373 ( <le-14) 
-0.141 (<le-14) 
-0.001 ( le+00 ) 
-0.342 (<le-14) 
-0.135 (<le-14) 
-0.040 ( 7e-02) 
-0.604 (<le-14) 
-0.383 (<le-14) 
-0.453 (<le-14) 
0.045 ( 4e-03 ) 
0.309 (<le-14) 
0.378 (<Ie- 14) 



0.708 (<le-14) 
0.745 (<le-14) 
0.720 ( <le-14) 
0.666 (<le- 14) 
0.396 (<le-14) 
0.585 ( <le-14) 
0.499 ( <le-14) 
0.321 (<le-14) 
-0.041 ( 6e-02 ) 
0.465 ( <le-14) 
0.164 (<le-14) 
0.107 (<le-14) 
0.701 (<le-14) 
0.435 (<le-14) 
0.619 (<le-14) 
-0.152 (<le-14) 
-0.288 (<le- 14) 
-0.553 (<Ie- 14) 



0.576 ( <le-14) 
0.530 ( <le-14) 
0.576 (<le- 14) 
0.557 (<le-14) 
0.334 (<le-14) 
0.386 (<le-14) 
0.362 (<le-14) 
0.307 (<le-14) 
-0.064 ( 3e-03 ) 
0.336 (<le-14) 
0.095 ( le-05 ) 
0.118 ( 6e-08 ) 
0.466 (<le-14) 
0.239 ( <le-14) 
0.477 ( <le-14) 
-0.163 (<le^l4) 
-0.117 ( 4e-14) 
-0.442 (<le-14) 



-0.073 ( 2e^06 ) 
-0.137 ( 3e-10) 
-0.115 ( 9e-14) 

0.003 ( 9e-01 ) 

0.047 ( 2e-03 ) 
-0.065 ( 3e-03 ) 
-0.065 ( 3e-05) 

0.101 ( 5e-ll ) 
-0.050 ( 2e-02) 

0.022 ( 3e-01 ) 
-0.016 ( 5e-01 ) 

0.064 ( 3e-03 ) 
-0.151 ( 3e-12) 
-0.088 ( 5e-05 ) 
-0.072 ( 3e-06) 
-0.075 ( le-06 ) 
0.165 ( <le-14) 
-0.080 ( 2e-04 ) 



" For each entry, we list the Spearman rank correlation coefficient (rg) and the probability of the null 
hypothesis ( Pnuii ) in parenthesis. 

* The black hole masses are calculated using the formalism presented in Vestergaard & Peterson 
(2006); Eddington ratios (L/LEdd) are calculated assuming that the bolometric luminosity Lboi ~ 
9ALa(5100). 

{X, Y;Z) denotes the partial correlation between X and F, controlling for Z. 



